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ABSTRACT 


Presented  are  the  results  of  studies  aimed  at  revealing  regularities  of  the  effect  of 
control  sensitivity  and  feel  system  characteristics  on  low-  and  high-frequency  PIO  and 
developing  methods  of  evaluating  this  eflfect.  The  study  has  been  performed  at  TsAGI 
in  compliance  with  the  AF  sponsored  research  contract  SPC-94-4027. 

The  technique  of  PIO  investigation  on  a  flight  simulator  has  been  refined.  Effect  of 
simulator  motion  on  PIO  has  been  studied.  It  has  been  shown  that  PIO  ground-based 
simulation  should  be  conducted  on  a  moving-base  simulator,  because  the  motion  cues 
can  influence  greatly  the  PIO  tendency  in  the  modeled  conditions  (in  some  cases  it 
becomes  greater,  in  others  it  decreases).  A  method  of  pilot  dynamic  performance 
identification,  which  allows  us  to  obtain,  from  the  experiments  directly,  a  common 
pilot  describing  function  and,  at  the  same  time,  dynamic  performance  of  closed  loop 
limb-manipulator  system. 

The  experiments  were  conducted  to  study  the  effect  of  feel  system  and  command 
sensitivity  characteristics  on  low-firequency  PIO.  Obtained  are  the  experimental  results 
concerning  the  optimum  feel  system  and  command  sensitivity  characteristics  for  a 
tracking  task,  and  pilot  ratings  worsening  in  the  case  of  the  characteristics  deviation 
from  their  optimum  values.  The  criterion  is  proposed  for  the  estimation  of  manipulator 
feel  system  and  command  sensitivity  influence  on  PIO.  The  criterion  allows  a  designer, 
firstly,  to  define  the  optimum  control  sensitivity  values  for  the  given  feel  system 
characteristics  for  a  tracking  task,  and,  secondly,  to  estimate,  in  terms  of  the  PIO 
tendency,  pilot  ratings  worsening  while  command  sensitivity  and  feet  system 
characteristics  deviating  from  their  optimum  values. 

High  frequency  PIO  occurring  for  certain  feel  system  characteristics,  aircraft  dynamic 
performance,  structural  elastic  mode  have  been  studied,  and  the  reasons  of  the 
phenomenon  arising  have  been  considered.  It  has  been  shown  that  for  an 
unmaneuverable  aircraft  one  of  the  causes  of  high  frequency  PIO  can  be  the  aircraft 
elastic  modes  coupled  with  limb-manipulator  dynamics.  The  degree  of  the  coupling 
depends  on  the  elastic  mode  amplitude  peak  vdue,  type  of  manipulator  and  its  feel 
system,  command  sensitivity  characteristics  and  aircraft  dynamic  performance.  The 
influence  of  accelerations  experienced  by  a  pilot  plays  a  significant  role  in  PIO 
phenomenon.  To  analyze  high-frequency  oscillations  in  the  roll  axis  the  mathematical 
pilot-aircraft  system  model  can,  with  a  good  reason,  be  used. 
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NOMENCLATURE 


F,X  -  manipulator  force  and  displacement  (kg,  mm) 

Fyi  ,  Xjj  -  longitudinal  command-response  gradients  (control  sensitivity  character¬ 
istics,  kg/g,  mm/g) 

Fp,Xp  -  lateral  command-response  gradients  (kg/deg/sec,  mm/deg/sec) 

F^ ,  -  feel  system  gradient  and  breakout  force  (kg/mm,  kg) 

n^,  Tiy  -  normal  and  lateral  accelerations 

% 

rir  -  normal  acceleration  per  unit  angle  of  attack  (g/rad) 

а,  p  -  angle  of  attack  and  sideslip  angle  (rad) 

б, <l>,y/  -  pitch,  roll  and  yaw  angles  (rad) 

<1  >P  >  f  ■  pitch,  roll,  yaw  rates  (rad/sec) 

V  -  flight  velocity  (m/sec) 

7^  -  roll  mode  time  constant  (sec) 

Yp,  -  transfer  functions  models  of  pilot  and  controlled  element 

Yfi  -  -  transfer  function 


2  ,Y  -  aerodynamic  forces  along  z-  and  y-  axes 

M  ,L,N  -  pitching,  rolling  and  yawing  aerodynamic  moments 

deflection  of  aileron,  elevator  and  rudder  (deg) 

7  Y  M  No,...  -  force  and  moment  aerodynamic  dimen- 

sionless  derivatives 
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INTRODUCTION 


Pilot-induced-oscillation  phenomenon  (PIO)  has  been  the  problem  of  considerable 
importance  for  modern  aircraft.  This  problem  crops  up  while  developing  almost  every 
new  aircraft.  However,  no  effective  design  or  experimental  methods  have  been  yet 
developed  to  predict  a  PIO  tendency  at  an  early  stage  of  aircraft  designing  and  tests. 
MIL-F-8785C  for  example,  just  states  that  an  aircraft  will  not  have  a  PIO  tendency, 
but  provides  no  guidance  in  the  area  of  precluding  PIO  tendency  by  design  So  the 
great  attention  has  been  paid  recently  to  investigation  of  this  problem  both  in  this 
country  and  abroad  . 

A  PIO  tendency  depends  on  aircraft  dynamic  performance,  as  well  as  on  manipulator 
feel  system  characteristics  and  command  -  response  gradients.  The  greater  success  was 
achieved  in  studying  the  effect  of  aircraft  dynamic  performance  (time  delay,  actuator 
rate  limiting,  etc.)  on  the  ordinary  low-frequency,  up  to  1  Hz,  PIO  (works  of 
D.McRuer,  T.Neal,  R.Smith,  R.Hoh^'^"'^  and  others).  For  studying  and  precluding  this 
type  of  PIO  there  are  well  developed  and  widely  used  methods  based  on  a  pilot-vehicle 
mathematical  model;  the  methods  allow  a  designer  to  estimate  the  dynamic 
performance  admissible  in  terms  of  pilot-vehicle  system  stability.  A  number  of  criteria 
are  developed  to  investigate  the  latent  PIO  causes  due  to  abrupt  changes  in  dynamic 
performance  or  in  flying  conditions  (presence  of  a  catalyst:  failures,  stress  situations, 
switching  from  one  control  loop  over  to  another,  etc.).  The  common  drawback  of  the 
approaches  is  that  they  do  not  take  into  account  the  effect  of  feel  system  and  control 
sensitivity  characteristics  on  PIO.  These  characteristics  are  usually  considered  to  be 
optimum  and,  as  a- rule,  they  are  considered  optimum  for  a  tracking  task.  However, 
neither  the  optimum  values  of  the  characteristics  nor  any  methodology  for  their 
definition  are  ever  shown.  These  drawbacks  diminish  considerably  the  approaches 
usefulness  since  in  reality  the  feel  system  characteristics  and  command  gradients  can 
greatly  differ  from  the  optimum  values. 

Along  with  low-frequency  PIO,  high-frequency  oscillations  (about  1.5-3  Hz),  named 
ratchet,  have  become  possible.  Especially  great  attention  has  been  paid  recently  to  the 
ratchet  in  the  roll  axis,  discovered  on  quite  a  number  of  aircraft  with  small  roll  mode 
time  constants.  As  it  is  shown  in  many  publications  (see  ^*’®^and  others),  ratchet,  as  well 
as  low-frequency  PIO,  occurs  as  a  result  of  pilot  /  aircraft  interaction.  The  main  cause 
of  ratchet  was  peaking  in  amplitude  of  a  limb-manipulator  describing  function  in  high 
frequencies  region.  The  presence  of  the  peaks,  their  values  and  frequencies,  depend  on 
manipulator  feel  system  and  command  sensitivity  characteristics.  However,  this 
dependence  and  a  ratchet  phenomenon  as  a  whole  have  not  been  sufficiently  studied 
yet.  Special  attention  should  be  paid  to  ratchet  coupling  with  structural  modes  of  an 
airframe  (CH-53,  F-111,  C-17  and  others),  since  there  is  a  tendency  in  modem 
aircraft  to  diminish  their  stmctural  stiffness  that  results  in  coupling  of  airframe 
stmctural  modes  and  limb-manipulator  dynamic  performance. 

A  number  of  publications  on  an  effect  of  manipulator  feel  system  and  command 
sensitivity  characteristics  on  handling  qualities  have  appeared  recently. 
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Nevertheless,  this  problem,  in  terms  of  PIO  tendency  especially,  is  insufficiently  studied 
yet,  and  at  present  there  are  no  methods  to  estimate  the  above  mentioned  effect. 

While  developing  controllability  criteria  and  theoretical  methods  of  PIO  tendency 
estimation,  it  is  also  important  to  improve  the  experimental  methods  of  investigation. 
Complexity  of  PIO  ground-based  investigation  is  determined  by  the  fact  that  aircraft 
oscillation  tendency  manifests  itself  irregularly  and  depends  not  only  on  aircraft 
characteristics,  but  on  the  piloting  task  and  pilot  physiological  state.  Therefore,  it  is 
necessary  to  improve  methods  of  PIO  modeling  and  experimental  data  processing  to 
study  both  evident  and  latent  causes  of  PIO. 

The  goals  of  the  work  are: 

•  development  of  PIO  investigation  methods  on  a  ground-based  simulator, 

•  studying  the  effect  of  manipulator  feel  system  and  command  sensitivity 
characteristics  on  low-ffequency  PIO  and  creating  the  criterion  of  the  effect 
estimation, 

•  studying  the  effect  of  limb-manipulator  d3mamic  performance  coupled  with 
structural  elastic  modes  on  high-frequency  oscillation  aircraft,  and  the  development 
of  a  technique  to  estimate  the  effect. 

Besides,  a  few  experiments  in  cooperation  with  MAI  and  FRI  were  conducted  on 
TsAGI  Flight  Simulator  FS-102  as  a  part  of  work  under  the  contracts  with  Wright 
Laboratory.  For  those  purposes,  the  Tu-154M  in-flight  simulator  dynamics  was 
modeled  in  FS-102,  the  tracking  task  indicator  was  modeled  on  the  special  "book-size" 
display,  the  flight, test  technique  was  worked  out  and  in-flight  simulator  model 
parameters  were  selected,  the  test-pilots  were  trained.  The  results  of  this  part  of  the 
work  are  not  considered  in  the  present  report,  since  they  are  referred  to  in  the  reports 
of  MAI  and  FRI. 
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Chapter  1 

EXPERIMENTAL  INVESTIGATION  TECHNIQUE 
1.1  Flight  Simulator 

The  experimental  part  of  the  work  was  conducted  on  flight  simulator  FS-102.  The 
simulator  is  intended,  mainly,  for  investigation  of  stability  and  controllability  of 
unmaneuverable  aircraft.  This  simulator  was  selected  because  its  structure  and  system 
parameters  give  an  opportunity  to  study,  more  completely,  the  effect  of  different  flight 
factors  on  PIO:  to  reproduce  linear  and  angular  accelerations  along  all  degrees  of 
freedom,  to  change  quickly  manipulator  types  and  their  characteristics  and  the  kinds  of 
flight  information  displayed,  and  other  flight  conditions. 

The  photos  of  the  simulator  and  its  systems  are  presented  in  fig.  1.1-1. 3.  The  simulator 
has  following  principal  characteristics: 

•  Visual  system;  single-channel,  optical  collimating  system,  computer-generated 
image  of  a  runway  and  its  vicinity  (fig.  1.2). 

•  Motion  system  (fig.  1 . 1):  of  synergetic  type,  6DOF  with  the  travel  limits; 

-  vertical  ±  1 .2  m,  longitudinal  and  lateral  directions  ±  1 .5  m; 

-  roll  ±30  deg,  pitch  ±  40  deg,  yaw  ±  60  deg. 

•  Pilot  cockpit  (fig.  1.2);  two  seats,  the  equipment  ordinary  for  unmaneuverable 
aircraft. 

•  Piloting  displays  (fig.  1.2);  In  the  instrument  desk  there  are  the  ordinary  piloting 
indicators  and  two  special  displays  installed.  On  the  first  of  them,  the  right  one, 
there  are  different  indicators  reproduced  which  display  the  current  values  of  angle- 
of-attack,  normal  acceleration,  airspeed,  altitude,  vertical  speed.  The  left  display 
(fig.  1.3)  was  used  as  an  indicator  for  the  tracking  task  simulation  (fig.  1.4)  while 
studying  PIO  phenomenon. 

•  Control  manipulators:  changeable.  The  spring  central  stick  and  electro-hydraulic 
side  stick  (fig.  1.3)  were  mostly  used  in  the  experiments  (the  characteristics  are 
presented  in  fig.  1.5,  1.6).  In  several  experiments  the  electro-hydraulic  central  stick 
with  widely  changeable  characteristics  was  used.  Using  these  manipulators  was 
determined  by  simultaneous  works  which  were  to  be  conducted  by  MAI  and  FRI 
using  these  types  of  manipulators. 


1.2  The  Technique  of  PIO  Simulation. 

1 .2. 1  Piloting  Task  and  Other  Experimental  Conditions. 

Piloting  task.  For  experimental  investigation  of  PIO  the  tracking  task  was  used  (the 
diagram  is  represented  in  fig.  1.4).  This  type  of  piloting  task  was  selected  due  to  its 
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methodological  advantages.  First,  only  in  the  case  of  persistent  handling  under 
disturbance  input  conditions,  pilot  describing  functions  can  be  identified.  It  is 
impossible  to  develop  theoretical  methods  of  controllability  or  investigate  PIO 
problems  without  these  functions.  Second,  it  was  this  very  piloting  task  that  was  used 
to  develop  a  number  of  existing  methods  to  study  controllability  as  a  whole  and  criteria 
for  PIO  tendency  evaluation  in  particular.  Therefore,  when  considering  this  piloting 
task,  there  is  a  possibility  to  compare  the  results  of  theoretical  and  experimental 
investigations;  it  results  in  improving  the  theoretical  methods  of  controllability 
investigation,  on  the  one  hand,  and,  on  the  other  hand,  in  simplifying  the  experimental 
results  analysis.  Third,  the  tracking  task  is  easy  to  simulate  and  it  is  easy  for  both  pilots 
and  operators  to  be  trained,  its  modeling  does  not  require  complex  experimental 
equipment  and  experimental  conditions;  results  obtained  in  different  works,  on  flight 
simulators  and  in  flying  conditions,  are  easy  to  reproduce  and  compare.  Finally,  the 
pilot-aircraft  system  responses  defined  for  this  piloting  task  can  be  a  sort  of  basic 
standard  for  PIO  tendency  evaluating.  This  very  task  was  used  to  conduct  unique  flight 
tests  within  different  programs  (LAHOS  and  others).  In  the  present  work  this  task  was 
used  as  well  for  the  joint  experiments  (TsAGI,  MAI,  FRI)  carried  out  to  develop  a 
technique  of  PIO  analyzing. 

In  addition  to  the  tracking  task,  a  landing  approach  was  modeled  in  the  course  of 
experiments.  The  experimental  conditions  for  a  landing  approach  task  correspond  to 
those  defined  in  the  report  and,  therefore,  in  this  report  are  not  referred  to. 

Trarkinp  task  indication  type.  While  performing  a  pitch  tracking  task,  a  pilot  was 
instructed  to  keep  the  pitch  tracking  error  e  (i.e.  the  difference  between  the  pitch 
predetermined  by  a.certain  disturbance  input  function  i  and  the  current  pitch  value  0  ) 
(fig.  1.4)  within  the  limits  represented  in  the  indicator  as  lines.  The  lines  were  to  be  the 
permissible  range  of  the  tracking  error.  An  indication  of  the  lines  and  the  tracking  error 
were  reproduced  in  the  attitude  indicator  together  with  an  indication  of  the  current 
values  of  pitch,  roll  and  yaw,  altitude,  vertical  speed  and  others.  Angles  of  pitch  and 
roll  were  reproduced  in  fiill  scale.  The  image  of  the  attitude  indicator  was  generated  on 
a  computer  and  then  reproduced  on  a  special  "book-size"  display  (fig.  1.3),  which  was 
installed  in  the  instrumentation  desk  of  the  simulator  cockpit.  Two  types  of  indication 
were  considered  (fig.l.7a,b).  In  the  first  case  (fig.  1.7a)  the  mark  of  the  tracking  error 
moved  while  the  lines  were  fixed.  In  the  second  case  (fig.  1.7b)  the  error  mark  was 
fixed  about  the  central  line  of  the  attitude  indicator  while  the  lines  moved,  i.e.  were 
deflected  in  accordance  with  the  pitch  angle. 

The  investigation  results  showed  that  pilots  adaptability  was  absolutely  the  same  in 
both  cases  of  error  indication.  Piloting  accuracy  does  not  depend  on  the  type  of 
indication  either.  However,  for  the  operators,  who  participated  in  the  experiments  as 
well,  the  first  type  of  indication  seemed  to  be  easier  to  comprehend  and  to  use  in 
practice.  Therefore,  the  first  indication  type  was  chosen  for  further  investigation 
(fig.  1.7a). 

To  simulate  a  roll  tracking  task,  the  roll  error  mark  was  presented  on  the  display  and 
the  error  was  to  be  nulled  by  a  pilot  while  performing  the  tracking.  For  this  case  no 
precision  limits  were  displayed. 
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Disturbance  input  and  aircraft  dynamics.  Pitch  and  roll  force  functions  manifested 
themselves  as  sum  of  sines  (SOS): 


i=15 

F(t)  =  X4  sin®,./ 

1=1 

Magnitudes  A-  and  frequencies  ta,  used  in  the  pitch  and  roll  loop  correspond  to  those 
plotted  in  fig.  1.8.  Higher  force  function  frequencies  for  the  roll  in  comparison  with  the 
pitch  are  accounted  for  by  the  fact  that  in  the  longitudinal  channel  low-frequency 
oscillations  were  studied,  while  in  the  lateral  channel  high-frequency  oscillations  were 
considered. 

In  a  few  experiments  some  LAHOS  dynamic  configurations  were  modeled  for  the 
results  to  be  later  compared  with  the  results  of  MAI  and  FRI  studies.  In  the  present 
work  LAHOS  configurations  1-4,  2-1,  2-10,  3-3,  4-10  and  others  were  considered. 

In  other  experiments  both  longitudinal  and  lateral  aircraft  motion  was  modeled  in 
accordance  with  the  equations  given  in  part  3.2. 

Three  test-pilots,  one  former  military  pilot  and  one  operator  participated  in  the 
experiments. 


1.2.2  An  Effect  of  Motion  Cues  on  PIO 


The  conducted  investigations  of  acceleration  effect  on  PIO  phenomenon  and  available 
publications  as  well,  show,  that  motion  cues  play  a  significant  role  in  PIO.  The  degree 
of  this  effect  depends  on  quite  a  number  of  factors:  control  channel,  aircraft 
performance,  flying  task  and  others. 


Let  us  consider  first  the  experimental  results  of  motion  cues  effect  on  roll  control,  see 
fig.l.9.-1.12.  Fig.1.9  illustrates  roll  damping  influence  (^Vj)  on  roll  tracking 

precision.  It  is  seen  that  simulator  motion  diminishes  a  roll  error  considerably.  The 
positive  role  of  motion  cueing  is  especially  evident  at  low  values  of  roll  damping, 
where  a  PIO  tendency  is  observed.  In  fig.  1.10  the  data  on  roll  rate  magnitudes  of  the 
oscillations  occurred  at  low  negative  roll  damping  values  are  shown.  The  data  were 
obtained  earlier  on  TsAGI's  simulator  with  and  without  motion  system  for  various 
fields  of  view  It  is  seen  that  angular  accelerations  influence  greatly  PIO.  Due  to  the 

simulator  motion  the  roll  rate  oscillation  magnitude  decreases  3-5  times. 


That  positive  effect  of  simulator  motion  on  handling  quality  can  be  illustrated  by  the 
pilot  describing  function  presented  in  fig.1.11.  The  data  show  that  in  the  case  of  both 
moving  and  unmoving  simulator,  the  pilot  behaviour  can  be  described  quite  well  by  the 
function 


Yp=Kp(TiS  +  l)e-^F 
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However,  as  a  simulator  moves  pilot  pure  time  delay  decreases  (in  this  case,  from 
r=  0.26  sec  to  r=  0.19  sec).  Due  to  time  delay  decreasing  a  pilot-aircraft  system 
stability  margin  increases,  that  allows  a  pilot  to  raise  his  gain  (in  our  case  from 
Kp-A.S  to  =  7.5).  For  this  reason  in  the  case  of  moving  simulator  roll  error 

compensation  improves  and  PIO  probability  decreases. 

In  some  cases  motion  cues  can,  on  the  contrary,  promote  PIO.  Roll  high  frequency 
oscillations  can  be  an  illustration  of  this  effect.  In  fig.  1.12  time  histories  obtained  on  a 
moving/unmoving  simulator  for  the  aircraft  model  with  certain  elastic  modes  are 
shown.  (The  results  of  the  experiments  will  be  considered  in  detail  in  chapter  3.)  It  is 
seen,  that  oscillations  appear  only  on  a  moving  simulator.  Neither  in  the  previous 
investigations  nor  in  the  present  one  the  attempts  to  reproduce  ratchet  cases  in  real 
flight  have  been  successful  on  unmoving  simulators.  This  fact  points  out  once  more  to 
the  essential  role  of  motion  cues  in  high  frequency  oscillations  phenomenon. 

In  the  longitudinal  channel  an  effect  of  motion  cueing  on  piloting  and  a  PIO  tendency  is 
not  so  evident.  Nevertheless,  according  to  the  comments  of  the  evaluation  pilots, 
simulator  motion  makes  simulation  conditions  seem  more  realistic.  The  pilots  noted 
that  in  some  cases  motion  cues  intensify  a  PIO  tendency,  in  others  they  do  not 
influence  the  tendency  or  mitigate  it.  Motion  cues  influence  piloting  precision  in  the 
same  way  as  they  affect  PIO. 

A  possibility  of  a  negative  effect  of  accelerations  on  PIO  is  consistent,  for  example, 
with  R.  Smith's  criterion  .  In  accordance  with  the  criterion  the  frequency  can  exist  at 
which  the  power  spectral  density  of  the  pilot's  normal  acceleration  due  to  pitch  attitude 
tracking  is  sufficiently  narrowband.  If  such  a  frequency  exists,  there  is  a  high 
probability  in  high  gain  tracking  task  pilot  will  switch  from  tracking  pitch  to  tracking 
the  normal  acceleration  he  feels  at  that  frequency.  If  phase  margin  of  the  pilot-felt 
normal  acceleration  to  stick  force  dynamics  is  less  than  zero,  then  the  aircraft  will  have 
a  tendency  to  PIO  at  that  frequency. 

The  positive  effect  of  cockpit  motion  on  piloting  is  usually  displayed  while  the 
accelerations  are  reproduced  with  regard  to  the  distance  between  the  cockpit  position 
and  center-of-gravity  1 »  0)  and  especially  while  there  is  no  PIO  tendency  observed. 

The  effect  of  cockpit  motion  on  pitch  tracking  may  disappear  at  1=0,  which  can  be 
attributed  to  the  fact  that  motion  cues  do  not  give  a  pilot  any  additional  information  in 
comparison  with  visual  ones.  It  has  been  shown  that  in  the  case  of  an  aircraft  of  a 
traditional  configuration  motion  cues  do  not  practically  lead  visual  cues  due  to  the 
inseparable  connection  of  a  normal  acceleration  and  a  pitch  angle  at  pilot  activity 
frequencies 


and  due  to  the  fact  that  the  acceleration  in  c.g.  does  not  significantly  lead  a  pitch  angle. 
Pitch  acceleration  feeling  and  its  using  by  a  pilot  are  hampered  due  to  the  strong  effect 
of  normal  accelerations  acting  in  combination  with  a  pitch  acceleration. 
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In  view  of  the  distance  between  the  pilot  cockpit  and  c.g.,  the 


-  transfer  function 


takes  the  form 


e 


n,  s 

‘-a 


I  . 
+— r . 
g 


It  is  seen  from  this  equation,  that  an  acceleration  leads  a  pitch  angle  and,  therefore, 
motion  cues  can  produce  a  favourable  effect  on  controlling.  If  the  aircraft  tends  to 
oscillate  in  pitch,  simulating  a  normal  acceleration  can  worsen  these  oscillations,  as  it 
was  observed  in  several  experiments. 

Thus,  the  data  considered  above  show  that  motion  cues  can  influence  piloting  and  a 
PIO  tendency  greatly.  Three  cases  should  be  defined:  motion  cues  producing  a  positive 
effect  on  piloting  and  mitigating  a  PIO  tendency;  motions  cues  producing  no  effect  at 
all;  motion  cues  intensifying  a  PIO  tendency.  At  present  there  are  no  theoretical 
methods  to  estimate  a  degree  of  this  effect  for  all  possible  cases.  Therefore, 
experimental  investigations  of  PIO  should  be  conducted  on  moving-base  flight 
simulators.  It  can  be  mentioned  also,  that  a  moving-base  simulator  has  some 
methodological  advantages  in  comparison  with  an  in-flight  simulator.  For  example,  on 
a  moving  -  base  simulator  there  is  a  possibility  to  change  acceleration  conditions 
(scaling,  separate  switching  on  different  degrees  of  freedom,  etc.)  without  changing 
other  flying  conditions,  which  is  impossible  in  real  flight  due  to  the  inseparable  unity  of 
motion  and  visual  cues  and  other  types  of  flying  information.  Therefore,  for  a  study  of 
motion  cues  effect' on  PIO  phenomenon,  a  moving-base  simulator  is  preferable  to  an 
in-flight  simulator. 


1.3  Experimental  Data  Processing 

In  the  work  both  subjective  and  objective  methods  of  experimental  data  processing 
were  used. 

Pilot  rating  scales  used.  In  the  work  two  pilot  rating  scales  were  used:  PIO  rating  scale 
(PIOR)  and  Cooper-Harper's  pilot  rating  scale  (PR).  The  PIOR  scale  was  used  to 
illustrate  that  pilot  ratings  worsen  due  to  PIO  tendency  intensification  as  command 
sensitivity  increases  or  feel  system  gradients  decrease  in  comparison  with  their 
optimum  values.  However,  for  the  final  rating  the  PR-scale  was  used,  since  it  is  a 
multipurpose  scale  adapted  for  evaluation  of  controllability  in  different  flight 
conditions,  including  a  PIO  tendency.  The  scale  is  well  known  for  experts  in  stability 
and  controllability  problems  and  pilots  engaged  in  the  experiments.  It  is  important  also 
that  this  scale  is  used  for  standardization  of  handling  qualities  in  Specifications  in 
different  countries. 

Each  studied  aircraft  configuration  was  flown  no  less  than  3  -  5  times.  The  ratings 
obtained  were  averaged.  In  accordance  with  the  technique  stated  in  the  confidence 
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interval  of  rating  arithmetic  mean  for  3-5  runs  does  not  exceed  PR  =  0.7-1.  For  the 
aircraft  configurations  of  Level  1,  the  confidence  interval  is  about  PR  =  0.5. 

The  technique  of  pilot  describing  function  identification.  In  the  study  the  run-time 
histories  processing  was  carried  out  according  to  the  technique  for  a  single-loop 
tracking  task  . 

In  this  study  an  attempt  was  made  to  develop  an  experimental  technique  of  an 
identification  of  a  pilot  describing  function  as  a  whole  and,  in  particular,  its  partial 
corresponding  to  the  dynamics  of  a  limb-manipulator  system.  The  block  diagram  of  the 
pilot-aircraft  system  for  this  case  is  shown  in  fig.  1.13.  The  model  does  not  contradict 
the  modern  pilot  behaviour  models,  for  example,  the  McRuer's  model  presented  in 
fig.  1 . 14.  In  this  diagram  and  further  the  following  notations  are  applied: 
e  -  tracking  error  displayed  on  an  indicator  and  observed  by  a  pilot  (command 
stimulus), 

i  -  disturbance  input  (sum-of-sines  forcing  function), 

6,(l>  -  current  state  variables  (for  example,  pitch  or  roll  angle), 

F  '  -  stick  force,  kg 

X  -  stick  displacement,  mm 

f  -  force  disturbance  generated  in  feel  system, 

-  transfer  function  model  of  central  nervous  system, 

7;^  -  transfer  function  model  of  closed-loop  limb-manipulator  system, 

-  transfer  function  model  of  neuromuscular  system, 

Yp  -  pilot  transfer  function  model  (J^  =  x  Yi^), 

Yfg  -  feel  system  transfer  function  model, 

-  transfer  function  model  of  the  controlled  element  (for  example,  stick 
displacement  {X,  mm)  referred  to  pitch  or  roll  (^,^ ,  deg)), 

-  pilot  remnant  transferred  to  visual  input, 

-  limb-manipulator  system  remnant  transferred  to  force. 

As  it  is  seen  from  the  diagram,  a  pilot-aircraft  system  incorporating  a  neuromuscular 
system  is  a  two-loop  model.  In  order  to  identify  simultaneously  two  transfer  functions 
in  thi.s  system  (Yp  and  its  partial  7/ni)>  the  inputs  i  and  /  should  be  uncorrelated.  Let 

us  consider  each  of  them  to  be  Gaussian  white  noise  passing  through  a  linear  filter. 

In  accordance  with  the  remnant  definition  given  in  remnants  and  are 

considered  uncorrelated  with  the  inputs  i  and  / . 

To  identify  the  transfer  function  models  the  Fourier  transform  algorithm  was  used. 
According  to  this  algorithm,  a  stick  displacement  and  a  tracking  error  can  be  described 
as  follows 

X{j(o)  =  X,  (ym)  +  {jco)  +  {jco)  +  {j(o) 

E(Jco)  =  {jco)  +  E„^  (Jco)  +  Ep  {jco)  +  E„^  (Jco) 

Subscripts  here  and  further,  refer  to  the  processes  in  a  closed-loop  pilot- 

aircraft  system,  caused  by  i{t),njt)j{t),nji). 
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With  the  use  of  the  transfer  functions  shown  in  diagram  1.13,  eq.(l.l)  takes  the 
following  form  (for  the  sake  of  brevity,  jco  is  omitted): 


Y  Y 

=  +  vv^ 

A  A  A(i +};,};,) 


F  + 


ns  js 


1  YY 

E  =  \l-^N.- 
A  A 


YfsY. 


(1.2) 


-F- 


(l  +  r„/^,)A  (1  +  44)  A 


N. 


where  A  =  1  +  4F^. 

Now,  taking  into  account  that  are  uncorrelated,  one  can  obtain  the 

following  power  spectra 


4 


9  - 

A(l  +  F^4) 


(1.3) 


where  8^1,3  -  power  spectra  of  inputs. 

Having  divided  left  ^d  right  hand  parts  of  the  first  two  equations,  one  has 


4(»  = 


40) 

40) 


(1.4) 


Taking  into  account  that 


S  =  F  +  f. 

Syi  , 

YfsUco)  =  -f- 


from  equations  (1.3),  it  is  easy  to  obtain 


40)  = 


444-4‘^.i4 

444 


(1.5) 


As  it  is  shown  in^'^\  for  a  two-loop  controlling  task  it  is  impossible  to  define  separately 
the  remnants  n^,n^  .  It  is  possible  to  define  their  sum  only.  Let  us  consider  that  the 
remnant  as  well  as  remnant  are  transferred  to  command  e .  So,  further  we  will 
consider  only  the  sum  of  remnant  ,  transferred  to  6 .  It  can  be  shown  that  the  power 
spectrum  of  the  remnant  can  be  calculated  as 
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(1.6) 


Se„e„{jo))  =  ■ 


SSl-SJl 


Thus,  in  the  case  of  a  two-loop  controlling  task,  the  describing  functions  and 

the  pilot  remnant  Se^e„  can  be  determined  by  (1.4),  (1.5),  (1.6) . 

As  to  the  practical  application  of  the  method,  the  following  fact  should  be  mentioned. 
A  pilot  describing  function  is  identified  more  precisely  when  an  input  function 
manifests  itself  as  sum  of  sines.  However,  in  the  case  of  a  two-loop  task,  the  using  of 
two  input  functions  leads  to  some  difficulties  in  the  identification  process.  These 
difficulties  are  accounted  for  by  the  difficulty  in  reproducing  two  uncorrelated  input 
functions.  Taking  this  into  account,  it  has  been  proposed  to  produce  one  of  the  inputs 
as  SOS,  and  another  input  as  white  noise  passing  through  a  linear  filter. 

The  first  experience  in  using  the  method  have  given  hopeful  results.  However,  the 
method  requires  further  developing. 
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Chapter  2 

EFFECT  OF  FEEL  SYSTEM  AND  COMMAND  SENSITIVITY 
CHARACTERISTICS  ON  LOW-FREQUENCY  P  I O 

Command  sensitivity  and  feel  system  characteristics  are  the  main  factors  affecting  PIO 
phenomenon.  This  fact  is  mentioned  in  a  number  of  publications.  It  is  enough  to  say, 
that  great  attention  is  paid  to  this  problem  in  Specifications,  for  example,  MIL.  It  is  for 
PIO  precluding  that  the  requirements  for  minimum  values  of  these  characteristics  are 
specified  in  them.  Nevertheless,  the  documents  available  show  the  effect  of  control 
sensitivity  and  feel  system  characteristics  on  PIO  insufficiently  and  in  kind  only.  It  is 
known  that  complex  interaction  of  such  factors  as  dynamic  performance,  feel  system 
characteristics  and  piloting  task  influences  a  degree  of  this  effect  considerably. 
However,  reliable  methods  of  the  effect  evaluation  have  not  been  developed  yet.  One 
of  possible  approaches  to  this  problem  is  considered  in  this  chapter. 


2.1  Effect  of  Command-Response  Gradients 

In  fig.2.1,  2.2  the  pilot  ratings  are  plotted  against  the  command-response  gradients  for 
longitudinal  and  lateral  channels.  The  relations  have  been  obtained  for  different 
dynamic  performance  for  an  aircraft  with  a  central  stick.  The  PIORs  presented  there 
together  with  the  PRs  show,  that  pilot  ratings  deterioration  in  the  case  of  command- 
response  gradients  'decreasing  is  correlated  with  PIO  tendency  intensified.  These  and 
quite  a  number  of  other  data  available  show  that  regularities  of  the  effect  of  control 
sensitivity  and  feel  system  characteristics  on  handling  qualities  and  PIO,  being  referred 
to  their  optimum  value,  are  the  same  for  different  piloting  conditions,  aircraft  classes, 
control  channels,  piloting  tasks,  dynamic  performance  and  manipulator  feel  systems 
(fig.2.3).  (The  above  mentioned  conditions  affect  only  optimum  values  of  command- 
response  gradients.  This  effect  will  be  considered  in  greater  detail  in  part  2.3.) 

As  it  is  seen  from  the  given  data,  an  increase  of  control  sensitivity  (command-response 
gradients  decrease)  over  optimum  values,  leads  to  PIO  tendency  arising.  At  small 
deviations  from  optimum  values  the  tendency  is  slight.  It  becomes  considerable  if 
control  sensitivity  characteristics  are  more  than  2  times  greater  in  comparison  with 
optimum  values.  For  example,  if  sensitivity  increase  is  2  times  the  PRs  and  the  PIORs 
deteriorate  by  0.6.  If  sensitivity  increase  is  4  times  the  PR-ratings  deteriorate  by  2-2.5, 
the  PIOR-ratings  deteriorate  by  2. 


dependence  of  pilot  ratings  on  dimensionless  command  gradients 

opt,...)  is  approximately  the  same  for  different  conditions  not  only  in 
P 

kind,  but  in  degree  as  well.  The  degree  of  pilot  ratings  deterioration  if  the  command- 
response  gradients  are  below  their  optimum  values,  can  be  described  by  the  relation 
(feel  system  characteristics  are  considered  optimum); 
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for  Cooper-Harper  scale 
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It  should  be  mentioned  also,  that  if  dynamic  performance  causes  no  PIO  tendency, 
control  sensitivity  decreasing  (command  gradient  values  increasing)  causes  only 
Cooper-Harper  pilot  ratings  variation,  while  PIO  ratings  do  not  change.  It  is  accounted 
for  by  the  fact,  that  command  gradients  increase  over  their  optimum  values  does  not 
result  in  PIO  tendency  intensification.  Controllability  worsening  in  this  case  is  due  to 
piloting  precision  deterioration  and  heavy  controlling,  but  not  to  PIO  tendency  (fig.2.1 
-  2.3).  If  dynamic  performance  is  a  cause  of  PIO  tendency,  control  sensitivity  reducing 
can  mitigate  it. 


It  is  evident,  that  eq.(2. 1)  describes  growth  of  PIO  tendency  only  approximately.  In 
each  particular  case,  PIO  severity  may  differ  from  this  empirical  relation,  as  it  is  shown 
in  fig. 2. 4.  This  tendency  could  be  evaluated  more  precisely  studying  pilot  behaviour 
models  in  terms  of  command  gradients  decrease.  Unfortunately,  the  models  and 
controllability  criteria  developed  so  far  (McRuer,  Neal-Smith,  et  al)  do  not  clearify  this 
tendency  pecularities.  In  fact,  they  assume  that  as  aircraft  gain  deviates  from  its 
optimum  value,  a  pilot  changes  his  gain  in  inverse  proportion  to  it  so,  that  YpY^ 

describing  function  remains  unchanged.  Therefore,  models  of  this  type  do  not  take  into 
account  the  fact  that  as  command  sensitivity  increases  a  pilot-aircraft  system  becomes 
more  oscillating.  To  define  possible  ways  of  eliminating  this  drawback  of  the  modern 
pilot  models,  special  investigations  of  pilot  adaptability  to  a  high  aircraft  gain  were 
conducted  in  the  present  work. 

In  fig.2.5  the  pilot  and  pilot-aircraft  describing  functions  obtained  for  the  tracking  task 
with  various  command  sensitivity  gradients  are  plotted.  The  analysis  of  these  and  other 
data  obtained  during  the  experiments,  shows  the  following:  if  control  sensitivity 
increase  exeeds  2  times,  a  pilot-aircraft  cut-off  (crossover)  frequency  becomes 
somewhat  higher,  i.e.  control  sensitivity  variation  interferes  with  pilot  adaptability.  The 
interference  becomes  considerable  at  high  frequencies.  As  it  is  seen  in  fig.2.5,  at 
frequencies  of  0,7-1  Hz,  the  magnitudes  of  pilot  describing  functions  are  practically  the 
same  in  spite  of  the  great  difference  in  aircraft  gains.  The  pilot  phase  remains  equal 
over  a  wide  frequency  range  for  all  aircraft  gains. 


This  pilot  peculiarity  may  be  a  result  of  different  muscles  of  an  arm  being  engaged  in 
deflecting  a  certain  manipulator  and  their  having  different  dynamics  and  displacement 
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ranges.  Wide-ranged  manipulator  displacements,  which  are  low-frequency  as  a  rule, 
are  produced  mainly  by  a  shoulder  and  a  forearm.  Due  to  this  fact,  pilot  adaptability  is 
higher  at  low  frequencies.  High-frequency  deflections  of  manipulator  are  produced  by 
muscles  of  a  hand,  having  more  narrow  displacements  limits.  It  is  these  limits  that 
restrain  the  pilot  high-frequency  adaptability  to  aircraft  gain  variations. 

Due  to  this  human  peculiarity  noticeable  changes  take  place  in  a  closed-loop  pilot- 
aircraft  describing  function.  As  it  is  seen  in  fig.2.5  ,  the  amplitude  peak  increases.  The 
peak  is  often  used  (Neal-Smith,  et  al)  as  a  standard  of  an  aircraft  tendency  to  PIO. 
Thus,  an  increase  of  control  sensitivity  leads  to  increasing  a  pilot-aircraft  cut-off 
frequency  and  rising  a  pilot  amplitude  ratio  in  a  high  frequency  band.  This,  in  turn, 
leads  to  peaking  in  a  pilot-aircraft  system  and,  therefore,  increasing  the  system 
oscillability.  This  is  a  possible  if  not  a  sole  cause  of  PIO  tendency  arising  under 
increased  control  sensitivity.  This  fact  should  be  paid  much  more  attention  to  in  the 
future. 


2.2  Effect  of  Feel-System  Characteristics 

The  main  feel  system  parameters  are  manipulator  gradient  F  and  breakout  force 
PIO  determined  by  abrupt  changes  in  force  gradient  values  have  been  considered  in  a 
number  of  works  (see  An  effect  of  force  gradients  and  breakout  on  PIO  is 

considered  in  the  present  work.  Its  regularities  have  not  been  sufficiently  presented  in 
publications  yet,  though  these  parameters  are  paid  great  attention  to. 

Fig.2.6;2.7  are  plots  of  pilot  ratings  and  force  gradients  for  different  force  breakouts. 
The  results  presented  were  obtained  for  both  central  and  side  sticks  for  optimum 
aircraft  control  sensitivity  characteristics.  These  and  other  data  available  brought  us  to 
the  following  results. 

Force  gradient  and  breakout  force  F*,  qualitative  influence  on  PIO  is  the  same  for 
different  manipulators,  aircraft  dynamics  and  piloting  tasks.  There  are  certain  optimum 
values  of  manipulator  force  gradients  for  each  manipulator  type  and  control  channel.  In 
the  case  of  any  deviations  of  force  gradients  from  their  optimum  values  the  aircraft 
tends  to  oscillate  and  pilot  ratings  worsen;  but  it  is  force  gradient  values  decreasing 
that  plays  a  decisive  role  in  PIO  tendency  severity  and  pilot  ratings  worsening.  For  low 
force  gradient  values  a  certain  additional  value  of  force  breakout  mitigates  a  PIO 
tendency  and  improves  controllability,  due  to  the  fact  that  too  low  gradient  values 
hamper  pilot's  measuring  control  forces.  In  this  case  the  lack  of  control  forces  is 
compensated  by  some  additional  breakout. 

It  should  be  mentioned  that  for  low  gradient  values  an  effect  of  control  sensitivity 
characteristics  on  PIO  is  also  greater,  what  has  been  shown  in  the  experiments,  see 
fig.2.8.  This  results  in  changing  relation  (2.1);  for  low  gradient  values,  deviations  of 
control  sensitivity  from  its  optimum  value  lead  to  more  considerable  pilot  ratings 
worsening  as  compared  to  (2. 1). 

If  gradients  and  breakouts  are  referred  to  their  optimum  values  (fig.2.9),  the 
dependence  of  pilot  ratings  on  these  referred  values  is  about  the  same  in  degree  for 
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different  manipulators  and  control  channels.  The  latter  fact  is  determined  by  the 
theoretical  approach  to  optimization  of  control  sensitivity  and  feel  system 
characteristics  presented  in  According  to  the  approach,  for  a  pilot  there  exist 
certain  desirable  ranges  of  force  F*  and  displacement  for  every  control 
manipulator.  The  optimum  gradient  value  (for  F^=0)  is  about 


(2.2) 


If  pilot  forces  and  displacements  differ  from  their  desirable  values,  controllability 
deterioration  degree  is  determined  by  and  •  This  dependence  remains  the 

same  for  different  manipulators  and  control  channels.  For  slight  deviations  of  and 

from  1  the  dependence  can  be  described  by  the  equation 


(2.3) 


where/  g  are  constant  for  different  manipulators. 

For  optimum  control  sensitivity  and  low  values  of  gradient  and  breakout  force,  a 
manipulator  displacement  range  does  not  depend  greatly  on  F  and  So,  it  can  be 
assumed  that 


•  (2.4) 

For  optimum  values  of  gradients  and  F^^=  0,  the  manipulator  force  range  is  close  to 
desirable  F^  .  If  gradient  values  decrease  and  breakout  values  increase  force  range 
values  change  according  to 


F  = 


•  F*  +  Fbr  . 


(2.5) 


Having  inserted  (2.5)  into  (2.3)  and  taking  into  consideration  (2.2)  and  (2.4),  we  have 


r 


APF  =  /Ig^ 


+ 


Fbr/ 

/  Fopt^*  j 


It  can  be  seen  from  the  equation  that  if  gradients  deviate  from  their  optimum  values, 
pilot  ratings  deterioration  is  determined  by  ^ and  ^y^pX-^  only,  regardless  of 

breakout  values  and  the  manipulator  type. 


Thus,  longitudinal  PIO  tendency  severity  due  to  deviations  of  feel  system  and  control 
sensitivity  characteristics  from  their  optimum  values  can  be  evaluated  with  the  function 
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plotted  in  fig.2.9.  Optimum  gradient  values  and  X*  for  different  manipulator  types  are 
presented  in  the  same  figure. 


2.3  Design  Criterion  for  Evaluation  of  Feel-System  and 
Command  Sensitivity  Characteristics  Effect  on  PIO 

An  essence  of  the  criterion  proposed  here  is  the  following; 

1.  With  the  help  of  the  criteria  available  (T.Neal,  R.  Smith,  R.Hoh  and  others)  severity 
of  PIO  tendency  is  evaluated  for  optimum  values  of  feel  system  and  control 
sensitivity  characteristics. 

2.  In  accordance  with  the  technique  stated  below,  the  tracking  task  optimum  values  of 
control  sensitivity  characteristics  are  evaluated  for  given  feel  system  characteristics 
and  aircraft  dynamic  performance. 

3.  With  the  help  of  the  function  in  fig.2.9  a  degree  of  PIO  tendency  increase  is 
evaluated  for  control  sensitivity  and  feel  system  characteristics  deviating  from  their 
optimum  values. 

To  develop  the  technique  of  optimum  values  selecting  for  control  sensitivity 
characteristics,  the  investigations  have  been  conducted  to  define  the  regularities  of  an 
effect  of  dynamic  performance  and  feel  system  characteristics  on  optimum  values  of 
command-response  gradients  for  the  tracking  task  conditions.  The  analysis  of  data 
obtained  shows  that  the  relation  of  optimum  command-response  gradients  to  aircraft 
dynamics  and  manipulator  feel  system  characteristics  is  the  same  in  its  nature  for  the 
tracking  task  and  for  the  ordinary  flying  task  (landing  approach,  cruise).  This  relation 
for  the  main  flying  tasks  was  considered  in  great  detail  in^^^^  and,  therefore,  in  this 
work  is  not  refferred  to.  To  illustrate  the  relation  of  optimum  command  gradients  to 
aircraft  dynamic  performance,  fig.2.10  presents  the  optimum  command-response 
gradients  which  have  been  obtained  in  the  experiments  for  five  different  LAHOS 

configurations.  Having  compared  configurations  2-1  and  5-11  which  have 
approximately  equal  damping  ratio,  one  could  see  that  if  natural  short-period  mode 
frequency  increases,  optimum  command-response  gradients  X^  increase  too.  The 

comparison  of  configurations  2-1  and  2-3  shows  that  damping  ratio  increase  results  in 
increasing  optimum  X  . 

The  experimental  data  obtained  allow  us  to  assume  that  optimum  command-response 
gradients  for  a  tracking  task  as  well  as  for  other  piloting  tasks  can  be  defined  by  A- 
criterion  .  For  the  longitudinal  channel  the  criterion  takes  the  form: 

1 

=  A(F^.Ft,r,..).  (2.6) 


1  + 


Vn 


2 

co^  + 


Tin 


\2 


V 


Yn  Uco*,xfi 
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To  show  that  A-criterion  application  is  well-grounded  for  a  tracking  task  as  well  as  for 
landing  approach,  aircraft  describing  functions  are  plotted  in  fig. 2.11.  These  functions 
were  defined  for  optimum  values  of  for  the  same  LAHOS  configurations  as 

presented  in  fig.2.10.  It  can  be  seen  that  in  spite  of  different  configurations  dynamics 
their  transfer  function  amplitude  ratio  curves  meet  at  about  the  same  point 
demonstrating  the  physical  nature  of  A-criterion  (2.6). 

But  quantitatively  optimum  values  of  command  gradients  depend  on  a  piloting  task. 
See,  for  example,  fig.2.10  where  the  optimum  command  gradients  for  some  LAHOS 
configurations  are  shown  for  the  tracking  task  and  landing  approach. 

It  obvious  that  the  difference  in  optimum  values  of  command-response  gradients  for 
different  piloting  tasks  can  be  considerable;  quantitatively  this  difference  depends  on 
aircraft  dynamic  performance. 

Relation  (2.6)  shows  optimum  command  gradients  values  in  kind  and  in  degree,  if  the 
parameters  in  (2.6)  are  specified.  The  comparison  of  the  calculated  and  experimental 

data  has  shown  that  the  values  of  the  parameters  in  (2.6),  but  for co^  do  not  depend  on 
a  piloting  task.  And  as  it  has  been  shown  in  ^  ^  the  value  of  To  can  be  put  equal  to 
140m/sec  for  all  flight  conditions  of  an  unmaneuverable  aircraft.  The  parameter  A 
depends  on  a  manipulator  type  and  its  feel  system  characteristics.  The  values  of  A  for 
central  and  side  sticks  of  an  unmaneuverable  aircraft  are  shown  in  fig.2.12. 

The  charateristic  frequency  (O*  depends  only  on  a  piloting  task  and  an  aircraft  class.  It 
is  1.2-1. 5  rad/sec  (or  0.7  rad/sec  for  a  landing  approach  task),  see  fig.2.11,  for  an 
unmaneuverable  aircraft  with  the  input  disturbance  shown  in  fig.  1.8. 

Thus,  optimum  command-response  gradients  values  for  a  tracking  task  on  an 
unmaneuverable  aircraft  with  central  or  side  sticks  can  be  calculated  using  (2.6),  where 

A  values  correspond  to  those  shown  in  fig.2. 12  and  co*  is  about  1 .2-1.5  rad/sec. 

Some  additional  studies  may  clearify  a  relation  of  parameter  A  to  manipulator  feel 

system  characteristics  and  specify  the  characteristic  frequency  co*  in  order  to  define 
optimum  command  gradient  values  for  other  control  channels,  manipulators  and 
aircraft  classes  for  a  tracking  task. 
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Chapter  3 


INVESTIGATION  OF  LIMB-MANIPULATOR  DYNAMIC  INTERACTION 
WITH  ROLL  CONTROL  OF  "ELASTIC”  AIRCRAFT 

3.1  Statement  of  The  Problem. 

In  view  of  change-over  to  fly-by-wire  system  on  modern  unmaneuverable  aircraft  a 
necessity  disappeared  for  a  pilot  to  apply  great  forces  to  move  a  control  linkage. 
Therefore,  small-inertia  manipulators  have  come  into  use  on  unmaneuverable  aircraft 
(small-mass  wheels,  mini-wheels,  central  and  side  sticks;  Tu-204,  A-320,  A-340,  etc.). 
For  these  types  of  manipulators  high  frequency  resonant  peaks  (1-3  Hz)  in  limb-ma¬ 
nipulator  system  describing  function  are  typical.  As  to  roll  mode  time  constants,  their 
low  values  are  not  characteristic  of  unmaneuverable  aircraft.  For  unmaneuverable 
aircraft  the  values  about  r^=0.5sec  or  more  are  typical.  Therefore,  for 
unmaneuverable  aircraft,  in  contrast  to  maneuverable  ones,  high  frequency  peaking  in 
limb-manipulator  system  and,  consequently,  ratchet  phenomenon  due  to  only 
unfavourable  limb-manipulator  system  characteristics  is  hardly  probable. 

But  for  a  modem  unmaneuverable  aircraft  noticeable  peaking  in  aircraft  describing 
function  in  a  frequency  band  about  2-3  Hz  is  typical  due  to  airframe  elasticity.  It  is 
accounted  for  by  the  fact  that  attempts  of  reducing  airframe  weight,  installing  engines 
on  a  wing,  increasing  aircraft  dimensions  lead  to  a  tendency  of  decreasing  aircraft 
elastic  mode  frequency  and  increasing  elastic  mode  amplitudes  on  modem  and 
prospective  aircraft.  As  a  result,  in  peak  frequencies  aircraft  dynamics  differs  from 
traditional  one,  which  could  be  described  by  roll  mode  only.  In  2-3  Hz  frequency  band 
a  resonant  peak  appears  in  roll  amplitude  ratio;  considerable  lateral  accelerations  can 
arise  while  roll  controlling.  (It  is  seen  from  flight  data  given  in  fig. 3.1  and  transfer 
function  models  responses  shown  in  fig.3.2a,b.)  Due  to  this  fact,  high-frequency 
oscillations  of  ratchet  type  become  possible  on  unmaneuverable  aircraft  as  well,  in 
spite  of  high  roll  mode  constant  values.  These  oscillations  were  observed,  for  example, 
in  flight  tests  of  one  of  Russian  unmaneuverable  aircraft  (fig.3.1). 

It  should  be  mentioned  that  the  phenomenon  in  question  may  present  a  problem  for  a 
maneuverable  aircraft,  since  high  frequency  oscillations  caused  by  stmctural  elasticity 
appear  on  maneuverable  aircraft  as  well.  This  type  of  oscillations  appeared,  for 
example,  on  a  relatively  light  military  aircraft  F-1 1 1  with  external  stores  loading  . 

This  chapter  of  the  report  aims  at  theoretical  and  experimental  validation  of  a 
possibility  of  roll  high-frequency  oscillations  caused  by  a  pilot  for  a  certain 
combination  of  manipulator  characteristics,  stmctural  elasticity  and  unmaneuverable 
aircraft  dynamics. 


21 


3.2  Unmaneuverable  Aircraft  Lateral  Motion  Model 
Coupled  with  Elastic  Modes 

Simplified  linearized  equations  of  aircraft  lateral  motion  coupled  with  structural  elastic 
modes  can  be  written  as  follows; 


<j)  - p-rXm6 

i+(D,+D,)i+(G  +  B,)4  =  R,S,  +  R,S,+li.S, 


where  ^  -  is  a  vector  of  aircraft  structural  modes  coordinates  of  dimension 

Nt,  where  Nt  is  a  number  of  structural  modes  considered, 
w  -  is  a  side  wind  gust  velocity. 

The  terms  describing  dynamic  interaction  of  structural  modes  and  rigid  aircraft  motion 
are  excluded  from  the  motion  equations  considered.  The  influence  of  static  elastic 
deformation  of  an  airframe  is  taken  into  account  by  means  of  special  corrections  of 
aerodynamic  coefficients  in  motion  equations  (3.1)  which  describe  an  aircraft  as  a  rigid 
body. 

The  matrix  of  aircraft  structural  stiffness  G  is  diagonal,  it  consists  of  squares  of 
structural  mode  natural  frequencies  in  vacuum;  the  matrix  of  structural  damping  is 
diagonal  as  well.  Matrices  G  and  do  not  depend  on  flight  conditions.  Matrices  of 
aerodynamic  stiffhess  B^,  aerodynamic  damping  and  "control  surfaces  efficiency" 
R^,R^,R^  change  due  to  flight  conditions  changing  and  depend  on  dynamic  pressure 
and  flight  velocity  Fin  a  first  approximation. 


Roll  rate  p  ,  yaw  rate  F  and  side  acceleration  affecting  a  pilot  or  registered  by 
control  system  sensors  (located  at  a  point  with  the  coordinates  xs,  zs)  are  computed 
according  to  the  following  equations; 

L  <Piix)^i( 


p{t,x)  =  p{t)  + 


lit) 


r[t,x)  =  r(t)  +  '^  \//i(x)^i(t) 


J=1 


X,  z)  =  — (t)  -  r(t)  cos  or  -  p(t)  sin  a]  -  «J(t) - 
--r(t)  +  -p(t)  +  -S[/,(^)  + 

g  g  g  ,=i 
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where 

fi(x) ,  <Pi{x) ,  y/^ix)  -  structural  i-mode  values  : 

fiix)  -  fuselage  deformation  along  the  y-axis,  m; 

(Piix)  -  fuselage  torsion  angle  with  respect  to  the  x-axis,  rad; 

y/i(x)  -  fuselage  bending  angle  with  respect  to  the  z-axis,  rad; 

xs,  zs  -  coordinates  of  the  point  considered,  m. 


A  first  approximation  of  elastic  aircraft  mathematical  model  takes  into  account  4 
structural  modes.  Matrices  for  cruise  conditions  (//=  11300m, 

M=0.S25,  ^1047kg/m^  ,  V-  243  m/sec)  for  aircraft  weight  83000kg  have  the  forms: 


1.49 

0.244 

-.178 

0.0 

47.2 

5.95 

-14.6 

0.0 

0.376 

0.212 

0.234 

0.0 

0.052 

3.73 

7.96 

0.0 

0.108 

0.306 

0.628 

0.0 

15.7 

5.24 

21.0 

0.0 

0.0 

0.0 

0.0 

0.26 

0.0 

0.0 

0.0 

6.26 

34.9 

69.5 

4555.7 

33.3 

rk 

1.20 

27.0 

63.5 

-28.1 

51.0 

29.4 

0.262 

37.0 

Diagonal  matrices  G,  are  the  following; 


173.4 

0.0 

0.0 

0.0 

0.211 

0.0 

0.0 

0.0 

0.0 

221.9 

0.0 

0.0 

,Dk  = 

0.0 

0.238 

0.0 

0.0 

0.0 

0.0 

244.3 

0.0 

0.0 

0.0 

0.250 

0.0 

0.0 

0.0 

0.0 

683.2 

0.0 

0.0 

0.0 

0.418 

As  a  result  for  cruise  flight  conditions  structural  modes  equations  can  be  written  as 
follows; 


IV 

1.703 

0.244 

-.178 

0.0 

'ir 

220.6 

5.95 

-14.6 

0.0 

[111 

1.2 

+ 

0.376 

0.450 

0.234 

0.0 

^2 

0.052 

225.6 

7.96 

0.0 

Is 

0.108 

0.0 

0.306 

0.0 

0.878 

0.0 

0.0 

0.678 

X 

^3 

+ 

-15.7 

0.0 

5.24 

0.0 

265.3 

0.0 

0.0 

689.5 

X 

I3 

LI4J 

LI4J 

34.9 

69.5 

45.7 

33.3 

1.20 

27.0 

63.5 

-28.1 

51.0 

29.4 

0.262 

37.0 

For  cruise  flight  condition  the  values  of  aerodynamic  derivatives  of  rigid  aircraft 
equations  considering  static  corrections  for  structural  elasticity  are  given  in  the 


lowing  table; 

derivative 

Y 

L 

N 

-0.1 

-3.75 

-2.2 

p 

0 

-0.9 

-0.076 

r 

0 

-0.3 

-0.33 

5a 

0 

-0.8 

0 

Sr 

-0.02 

-0.75 

-1.05 
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Values  of  structural  modes  in  different  fuselage  points  are  equal  to 

•  in  a  cockpit 

/  =  [-0.155  0.219  -0.293  -0.97 1] 

^  =  [0.0374  0.0148  -0.021  0.0027] 

=  [0.0093  -0.0092  0.0276  0.109] 

•  in  a  control  system  sensors  location 
/  =  [-0.0323  0.116  0.035  0J35] 

(^  =  [0.0326  0.0107  -0.0185  0.00455] 

=  [0.0048  -0.00015  0.016  0.017] 


transfer  functions  models  responses  measured  in  a  cockpit  are 
presented  in  fig.3.2.a,b. 


Block-diagrams  of  lateral  control  system  selected  for  simulation  are  given  in  fig.3.3. 


Feed-back  and  feed-forward  gains  of  control  system  for  cruise  flight  conditions  are 
equal  to  the  following; 


=  3.0  sec;  =  0.087  .^:^=0.6sec;  K,  =  1.0  sec; 

K.  -  0.0017  K,  =  0.0052  K,  =  0.003 


Coordinates  of  control  system  sensors  are  xs  =  4.6  m,  zs  =  0. 

Longitudinal  aircraft  motion  model  for  simulation  is  described  with  the  following 
equations; 


d  =  -Z^cc  +  q  -  Zg  Sg 
J  e 

q  =  M(^a  +  Mqq  +  Mgd  +  Mg  Sg 

Vi  1  X 

-  d{t)]  + 


Aerodynamic  coefficients  accounting  for  static  corrections  for  structural  elasticity  are 
given  in  the  following  table; 


Transfer  functions  models  of  side  acceleration  and  roll  rate  referenced  to  an  aileron 
deflection  take  the  forms: 

riy  (5+ 1.62)(/ +0,885  +  1.59)  (5- 7.2)(5  +  5.5) 

—  =  -0.187  (^  +  o.83)(5^ +0.545  +  2.2)  (5^  +  1.65  +  219)^ 

(5^+0.715  +  227)  (5^-6.35  +  555)  (5- +4.45  +  631) 

^(5^+0.415  +  222)  (5^+5+270)  (5“ +  0.635  +  678) 
p  2.44  (5^+0.445+2.27)  (5- +0,215+221)  (5- + 1.195  +  250) 

(5  +  0.83)  (5- +0.545  +  2.2)  (5' + 1.65  +  219)  (5' +0,415+222)  ^ 

(s  +  9.3)(s-8.66)  (5^+0.6785  +  690) 

^  (5^+5  +  270)  (5^+0.635  +  678) 


3.3  Experimental  Results 

The  results  of  high-frequency  pilot  assisted  oscillations  modeling  are  considered  in  this 
part  of  the  work. 

The  experiments  modeling  aircraft  dynamics  on  a  flight  simulator  (see  part  3.2)  show 
the  following: 

1.  Aircraft  elastic  mode  influence  high-frequency  oscillations  in  pilot-aircraft  system  to 
a  considerable  extent. 

High-frequency  oscillations  were  regularly  observed  in  the  course  of  experiments  on 
the  moving-base  simulator  for  both  step  roll  manipulator  input  and  permanent 
manipulator  deflections  determined  by  a  tracking  task.  This  fact  was  mentioned  by  all 
the  three  pilots  who  took  part  in  the  experiments.  In  fig.3.4,  3.5  the  time  histories  for 
an  elastic  and  rigid  aircraft  are  given;  the  moving-base  simulator  was  flown  by  one  of 
the  pilots.  The  comparison  of  the  given  time  histories  shows  that  high-frequency 
oscillations  appeared  on  the  elastic  aircraft  only. 

The  oscillations  frequency  remained  about  2.5  Hz  in  all  the  cases  regardless  of  a  pilot 
and  his  piloting  manner.  This  frequency  coincides  with  that  obtained  during  flight  tests 
(fig.3.1)  and  corresponds  to  the  frequency  of  the  first  aircraft  elastic  mode  (fig. 3. 2a, b). 
The  latter  fact  proves  it  was  aircraft  elasticity  that  caused  oscillations  in  both  ground- 
based  and  in-flight  experiments. 

2.  The  magnitude  of  such  high-frequency  oscillations  depends  on  a  resonant  peak 
magnitude  in  elastic  aircraft  transfer  function  responses. 

It  can  be  seen  in  fig.3.5-3.8  that  as  the  first  aircraft  elastic  mode  peak  becomes  higher, 
the  magnitudes  of  lateral  and  roll  rate  oscillations  increase,  the  magnitudes  being 
proportional  to  the  elastic  mode  peak,  see  fig.3.9. 
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As  it  is  shown  in  fig. 3. 10  the  pilot  ratings  worsen  as  the  peak  magnitude  increases.  The 
threshold  peak  magnitude  which  corresponded  to  ratchet  arising,  was  about  12,5%  of 
the  peak  maximum  value,  according  to  the  pilots.  As  the  peak  became  higher,  the 
pilots  noticed  the  unfavourable  effect  of  the  accelerations.  The  attempts  to  counteract 
these  disturbances  failed.  Smooth  manipulator  resetting  to  the  neutral  position  or 
setting  it  free  damped  these  oscillations. 

3.  Pilot-felt  accelerations  influence  high-frequency  oscillations  considerably. 

This  follows  from  the  comparison  of  the  data  obtained  in  the  experiments  while 
different  degrees  of  freedom  were  engaged:  roll  and  lateral  displacement  (fig. 3. 5),  roll 
only  (fig.3. 1 1),  no  degrees  of  freedom  switched  on  (fig.3. 12). 

The  data  showed  and  the  pilots  noticed  that  aircraft  oscillations  arose  on  a  moving 
simulator  only.  An  aircraft  oscillation  tendency  was  observed  even  while  only  the  roll 
degree  was  switched  on,  but  it  was  lateral  accelerations  that  influenced  this  tendency 
greatly.  An  acceleration  effect  interfered  with  tracking  task  performing,  thus,  the 
piloting  precision  worsened. 

4.  The  high-frequency  oscillations  are  characteristic  of  both  central  and  side  sticks  in 
spite  of  the  latter  having  an  armrest  and  a  damping  device,  see  time  histories  in 
fig.3.5-3.7  for  the  central  stick  and  fig.  3.13,  3.15  for  the  side  stick. 

A  tendency  to  high-frequency  oscillations  depends  on  the  direction  of  forces  applied. 
This  scewness  depends  on  the  manipulator  type.  In  the  case  of  a  central  stick  the 
tendency  is  greater  .while  the  stick  is  deflected  to  the  right,  but  for  a  side  stick  it  is  vice 
versa,  see  fig.3. 5  and  3.13.  According  to  the  pilots,  difference  in  the  type  of  scewness 
due  to  the  type  of  a  stick  is  accounted  for  by  the  fact  that  different  muscle  groups  are 
engaged  in  controlling,  these  muscles  having  different  dynamic  and  force 
characteristics.  In  the  case  of  a  central  stick  an  arm  and  upper  body  are  engaged,  while 
in  the  case  of  a  side  stick  with  an  arm  on  an  armrest,  only  a  forearm  and  a  hand  are 
used. 

5.  Aircraft  command  sensitivity  affects  high-frequency  oscillations  caused  by  structural 
elasticity  to  a  considerable  extent,  which  can  be  seen  from  the  time  histories  (fig.3. 5, 
3.14,  3.15)  and  from  the  relation  of  pilot  rating  to  command  sensitivity  as  well 
(fig.3.16). 

It  should  be  mentioned  that  pilot  ratings  variation  is  accounted  for  by  ratchet  mainly, 
according  to  the  pilots;  if  there  is  no  ratchet  observed  aircraft  gain  variation  does  not 
influence  pilot  ratings.  According  to  our  data  if  aircraft  gain  is  2  times  less  in 
comparison  with  its  optimum  value,  there  are  practically  no  high-frequency  oscillations 
observed. 

This  peculiarity  of  command  sensitivity  effect  on  high-frequency  oscillations  caused  by 
structural  elasticity  is  in  agreement  with  the  data  given  in  concerning  a  command 
sensitivity  effect  on  ratchet  at  low  roll  mode  time  constant  values. 
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6.  High-frequency  oscillations  are  possible  to  simulate  on  a  ground-based  simulator 
with  a  motion  system.  It  follows  from  all  the  data  presented  above. 

According  to  the  pilot  who  took  part  in  the  in-flight  experiment  (see  time  histories  in 
fig. 3.1)  and  in  the  ground-based  experiments  (see  time  histories  in  fig. 3. 17)  high- 
frequency  oscillations  felt  on  the  ground-based  and  in-flight  simulators  are  basically  the 
same.  However,  the  oscillation  tendency  observed  in  flight  experiments  was  less.  It  is 
accounted  for  by  two  facts:  first,  a  miniwheel  was  used  in  flight  while  central  and  side 
sticks  were  used  on  a  simulator;  second,  the  command  sensitivity  characteristics 
differed. 

It  has  been  mentioned  in  some  works  (see,  for  example,  ^’^)  that  ratchet  caused  by  low 
roll  mode  time  constant  values  is  difficult  to  reproduce  on  a  ground-based  simulator. 
As  to  high-frequency  oscillations  caused  by  structural  elasticity,  our  experience  shows 
that  this  type  of  oscillations  is  easily  reproduced  and,  thus,  can  be  studied  on  a 
simulator. 

The  experimental  data  concerning  pilot  describing  functions  will  be  considered  in  the 
next  part. 


3.4  Analysis  of  Pilot- Aircraft  System  Characteristics 

To  reveal  causes  of  high-frequency  PIO  and  ways  to  preclude  this  phenomenon  let  us 
consider  the  describing  functions  of  a  pilot  and  an  open-loop  pilot-aircraft  system 
obtained  in  the  co.urse  of  the  experiments,  see  fig.3. 18-3.21.  In  these  figures  the 
approximation  of  these  describing  functions  in  terms  of  a  pilot  transfer  function  model 
is  shown  as  well.  The  pilot  transfer  function  is  as  follows 

=  (3-2) 

where 

Kj,,Tj  -  pilot's  gain  and  lead, 

T  -  pilot's  equivalent  time  delay  (as  a  combination  of  pilot's  pure  time  delay 
and  computer  time  delay). 

In  accordance  with  the  diagram  in  fig.  1.1 4  -  transfer  function  model 

(7^  =  0,  =  00  )  can  be  presented  as 


+  2^1  T,s  +  \T^s'^  +2^ +  \  ^  ’ 

Let  us  consider  the  pilot/pilot-aircraft  describing  functions  in  fig. 3. 18  and  3.19  for  the 
case  of  roll  motion  described  with 

(3-4) 
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It  can  be  seen  that  there  is  noticeable  peaking  in  amplitude  ratio  at  high  frequencies  (1- 
3  Hz)  which  is  determined  by  limb-manipulator  system  dynamics.  As  a  result  a 
resonant  peak  in  an  open-loop  pilot-aircraft  system  appears  at  these  frequencies  if 
values  of  roll  mode  time  constant  are  low  {Tr  =  0.1  sec,  fig.3.18).  The  same  type  of 
peaking  was  observed  by  other  investigators  who  analyzed  ratchet  occurring  on 
maneuverable  aircraft  in  real  flight.  Some  studies  have  shown  (see  for  example) 
that  this  peak  magnitude  is  a  measure  of  high-frequency  PIO  tendency  caused  by  low 
roll  mode  time  constant  values.  Roll  high-frequency  oscillations  arise  when  a  peak 
magnitude  is  about  -6db  or  more.  Thus,  the  more  the  peak  magnitude,  the  greater  the 
PIO  tendency. 

If  roll  mode  time  constant  exceeds  0.5  sec  which  is  typical  of  “rigid”  unmaneuverable 
aircraft,  resonant  peaks  are  below  -6db.  This  can  be  seen  from  the  describing  functions 
presented  in  fig.3.19  and  3.20;  in  the  first  case  roll  motion  corresponded  to  eq.(3.4) 
where  Tr  =  0.5  sec;  in  the  second  case  roll  motion  corresponded  to  eq.  (3.1)  where 
Tr  =  1.2  sec  (structural  elasticity  was  not  taken  into  account).  Thus,  for  a  rigid 
unmaneuverable  aircraft  no  ratchet  was  observed  neither  in  our  experiments  no  in 
other  studies. 

The  experiments  were  conducted  for  an  elastic  unmaneuverable  aircraft  with  various 
elastic  mode  amplitudes  and  control  sensitivity  characteristics.  The  experiments 

showed  that  if  a  resonant  peak  in  exceeded  -6db,  high-frequency  PIO  arose.  The 

pilot/pilot-aircraft  describing  functions  are  given  in  fig.3.22  for  the  case  when 
structural  elastic  modes  (see  eq.(3.1))  corresponded  to  real  ones.  (As  it  has  been 
mentioned  above,  in  this  case  the  PIO  tendency  was  extreme.)  It  is  seen  that  for  this 
case  the  resonant  peak  magnitude  was  up  to  +6db. 

These  data  has  shown  that  the  pilot  and  his  neuromuscular  system  describing  functions 
obtained  in  the  experiments  can  be  adequately  approximated  by  transfer  functions 
models  (3.2),  (3.3). 

The  above  mentioned  facts  allow  us  to  conclude  that  both  high-frequency  oscillations 
on  unmaneuverable  aircraft  due  to  their  structural  elasticity  and  ratchet  phenomenon 
on  maneuverable  aircraft  due  to  low  values  of  roll  mode  time  constant  can  be  studied 
on  the  basis  of  pilot-aircraft  model  responses. 

The  pilot-aircraft  model  responses  can  be  defined  empirically  or  derived  from  pilot  and 
limb-manipulator  transfer  functions  models  (3.2),  (3.3).  In  the  latter  case  a  pilot's  pure 
time  delay  can  be  assumed  constant,  for  example,  t=0.3  sec;  parameters  Kp  and  Tj  can 
be  derived  from  "crossover  model"  as  the  data  in  fig.3.18-3.21  prove: 

^  JO) 

where  cOc  -  pilot-aircraft  crossover  frequency  (TpF^Ocu^)  =  1)  which  does  not  depend 
on  aircraft  characteristics. 
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It  should  be  mentioned  that  a  limb-manipulator  model  has  not  been  sufficiently 
developed  yet.  There  is  no  clear  idea  of  the  order  of  a  limb-manipulator  model  which 
would  be  adequate  to  analyze  different  cases  of  high-frequency  PIO.  It  was  concluded 
in  that  a  third-order  transfer  function  model  is  enough  to  describe  ratchet  caused  by 
low  roll  mode  time  constant.  The  pilot  describing  functions  presented  in  fig.3. 18-3.21 
show  that  there  are  two  resonant  peaks  in  pilot's  amplitude  ratio  at  the  frequencies 
exceeding  1  Hz.  These  describing  functions  can  be  matched  to  a  transfer  function 
model  of  the  fourth  order.  Comparing  the  experimental  and  calculated  data  showed 
that  an  adequate  approximation  of  pilot-aircraft  describing  functions  as  well  as  the 
resonant  peaks  is  achieved  if  transfer  function  model  (3.3)  is  applied.  We  may  conclude 
that  while  studying  high-frequency  PIO  a  limb-manipulator  transfer  function  such  as 
(3.3)  can  be  with  a  good  reason  applied. 

It  has  been  mentioned  in  the  publications  that  parameters  depend  on  a 

manipulator  type  and  its  feel  system  characteristics.  They  can  be 
7^2  =  0.05  -  0.15  sec,  ^  =  ^.05  -  1.0. 

In  our  experiments  these  parameter  magnitudes  were  Tj  =  0.12  sec 
=  0.2  sec,  T2  =  0.055  sec,  ^2  =  0.1.  However,  these  parameters  adjustment 
rules  determined  by  a  manipulator  type  and  feel  system  characteristics  have  not  been 
developed  yet.  This  hampers  an  application  of  mathematical  models  of  a  pilot-aircraft 
system  for  high-frequency  PIO  analysis.  These  rules  are  also  necessary  for  describing 
limb-manipulator  system  dynamics  while  developing  controllability  criteria  and 
analyzing  low-frequency  PIO.  Further  studies  to  develop  such  adjustment  rules  should 
be  carried  out. 

The  pilot  mathematical  model  considered  here  disregards  pilot-felt  lateral 
accelerations.  Therefore,  influence  of  these  accelerations  on  low-  and  high-frequency 
PIO  should  be  studied  in  greater  detail. 
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CONCLUSIONS 


1.  The  technique  of  PIO  experimental  studies  has  been  improved  while  modeling 
compensatory  tracking  task  on  a  ground-based  simulator. 

It  has  been  proved  that  ground-based  PIO  studies  should  be  conducted  on  a 
moving-base  simulator,  since  motion  cues  felt  by  a  pilot  can  influence  a  PIO 
tendency  greatly  (in  some  cases  the  tendency  is  mitigated,  in  others  it  is 
intensified). 

The  technique  has  been  improved  to  define,  from  experiments  directly,  a  common 
pilot  describing  function  and,  at  the  same  time,  dynamic  performance  of  closed 
loop  limb-manipulator  system.  The  technique  suggests  the  use  of  two  uncorrelated 
disturbance  functions:  one  of  them  is  the  visual  disturbance  function,  and  the  other 
is  the  manipulator  force  disturbance  function  generated  in  a  feel  system. 

2.  Main  regularities  of  a  feel  system  and  control  sensitivity  characteristics  effect  on 
PIO  have  been  revealed. 

It  has  been  shown  that  a  decrease  of  command-response  gradients  and  force 
gradient  in  comparison  with  their  optimum  values  leads  to  PIO  tendency 
intensification.  Greater  values  of  breakout  forces  at  low  values  of  force  gradients 
mitigate  PIO  tendency. 

As  control  sensitivity  increases  pilot  gain  adaptability  to  aircraft  gain  variation  is 
upset.  First,  the  pilot-aircraft  model  cut-off  frequency  becomes  somewhat  higher. 
Second,  at  the  frequencies  exceeding  the  cut-off  value,  pilot  model  amplitude  ratio 
curves  tend  to  converge.  The  amplitudes  are  practically  the  same  when  the 
frequency  is  about  0.7-1  Hz,  in  spite  of  the  great  difference  in  aircraft  gains.  The 
pilot  phase  remains  about  the  same  for  all  aircraft  gains.  This  is  accounted  for  by  the 
fact  that  quite  a  number  of  muscle  groups  with  various  dynamic  responses  and 
displacement  limits  participate  in  deflecting  a  particular  manipulator.  This  pilot 
peculiarity  should  be  studied  in  greater  detail  further,  as  it  could  form  the  basis  for  a 
mathematical  method  to  evaluate  the  control  sensitivity  effect  on  PIO. 

The  criterion  is  proposed  for  an  estimation  of  manipulator  feel  system  and  command 
sensitivity  influence  on  PIO.  The  criterion  allows  a  designer,  firstly,  to  define 
optimum  control  sensitivity  values  for  a  tracking  task  for  the  given  feel  system 
characteristics,  and,  secondly,  to  estimate,  in  terms  of  PIO  tendency,  pilot  ratings 
worsening  in  the  case  of  characteristics  deviation  from  their  optimum  values. 

3.  It  has  been  shown  that  there  is  a  tendency  to  high  frequency  oscillations  on  an 
unmaneuverable  aircraft  with  certain  characteristics  of  a  manipulator,  structural 
elasticity  and  aircraft  dynamics. 

High  frequency  PIO  is  possible  to  imitate  on  a  moving-base  simulator.  According  to 
the  pilot  the  high  frequency  oscillations  felt  on  ground-based  and  in-flight  simulators 
are  essentially  the  same. 
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Aircraft  elastic  modes  affect  considerably  high  frequency  oscillations  in  a  pilot- 
aircraft  system.  Their  mere  arising  and  their  intensity  depends  on  the  magnitude  of 
the  first  aircraft  elastic  mode.  Their  frequency  corresponds  to  the  frequency  of  the 
first  elastic  mode. 

High  frequency  PIO  intensity  depends  on  limb-manipulator  system  characteristics  as 
well;  these  characteristics  are  determined  by  the  manipulator  type  and  its  feel  system 
parameters.  High  frequency  PIO  appeared  in  the  experiments  with  both  central  and 
side  control  sticks.  High  frequency  PIO  tendency  depends  on  the  direction  of  forces 
applied.  This  skewness,  in  its  turn,  depends  on  the  manipulator  type.  In  the  case  of  a 
central  stick  the  tendency  is  greater  in  the  right  direction,  in  the  case  of  a  side  stick 
it  is  greater  in  the  left  direction. 

Aircraft  command-response  gradients  influence  high  frequency  oscillations  caused 
by  structural  elasticity  to  a  considerable  extent.  As  aircraft  gain  increases  and  pilot 
adaptability  gets  upset,  a  resonant  peak  in  an  open  loop  pilot-aircraft  system  is 
higher  and,  therefore,  a  PIO  tendency  is  intensified. 

Roll  and  lateral  accelerations  felt  by  a  pilot  play  an  important  role  in  high  frequency 
oscillations  phenomenon.  High  frequency  oscillations  were  observed  only  on  the 
simulator  with  a  moving  base. 

Causes  of  high  frequency  PIO  and  ways  of  its  precluding  can  be  adequately  studied 
by  means  of  analyzing  pilot-aircraft  model  describing  functions. 
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Fig. 1.5.  Central  Stick  Feel  System  Characteristics 


Fig. 1.6.  Side  Stick  Feel  System  Characteristics 
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Fig.  1.8  Forcing  function  spectral 

densities. 
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.  Lateral  Comiand-Response  Gradients  Influence  on  Pilot  Ratings 
(PR  and  PIGR)  for  two  values  of  roll  lode  time  constant  for 
an  aircraft  with  a  central  stick. 
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Fig. 2, 9  Pilot  Ratings  Deterioration  Degree  due  to  PIO  Tendency 
Increasing  in  the  case  of  Feel  Systei  (a)  and  Control 
Sensitivity  (b)  Characteristics  Deviation  froi  Their 

Optiittui  Ualues. 
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Fig. 2. 10  Optilui  Coiniand-Response  Gradients  for  5  LAHOS 
Configurations  for  Landing  Approach  and  Tracking  Task. 
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Fig.2.11  Roll  transfer  function  responses  for  optimal  control 
sensitivity  of  different  dynamic  configuration, 
al  tracking  task 
b)  landing  approach. 
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Fig.2*i2  Paraieter  A  Depending  on  Gradient  of 
Central  and  Side  Sticks  (Fbr=Ffr-F^  -0) 
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Fig.  3.1.  Time  HLsForLe^  foF  ihe  l^a-tcFei  Case. 
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Fig.3.2aRoll-rate  transfer  function  responses  for  the  aircraft 
with  .and  without  airframe  easticity. 


Fig. 3, 9*  Elastic  Mode  Aiplitude  A  Influence  on  Roll  Rate  Ap 
and  Lateral  Acceleration  Aj^,^  Amplitudes. 
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Histories  of  Elastic  Aircraft 


Manipulator  *  central  sticl:: 
Aircraft  gain:  K/Kn  =  1-5 


Elastic  mode  amplitude:  A  :::  1 
Motion  simulated:  roll 


ilot/Pilot— Aircraft  describing  function 
(isolated  roll,  Tr  =■  0.1s) 


Fig.3.20.Pilot/P] lot— Aircraft  describing  function 
(isolated  roll,  Tr  =  ,1,0s) 


